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Abstract 

This paper describes the applications of the fractured continuum model to the different enhanced geothermal 
systems reservoir conditions. The capability of the fractured continuum model to generate fracture characteristics 
expected in enhanced geothermal systems reservoir environments are demonstrated for single and multiple sets of 
fractures. Fracture characteristics are defined by fracture strike, dip, spacing, and aperture. The paper demonstrates 
how the fractured continuum model can be extended to represent continuous fractured features, such as long 
fractures, and the conditions in which the fracture density varies within the different depth intervals. Simulations of 
heat transport using different fracture settings were compared with regard to their heat extraction effectiveness. 
The best heat extraction was obtained in the case when fractures were horizontal. A conventional heat extraction 
scheme with vertical wells was compared to an alternative scheme with horizontal wells. The heat extraction with 
the horizontal wells was significantly better than with the vertical wells when the injector was at the bottom. 

Keywords: Geothermal reservoir simulation; Enhanced geothermal systems; Heat extraction: Fracture network; 
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Introduction 

The major objective of this work is to demonstrate the ap- 
plications of the fractured continuum model (FCM) to 
common enhanced geothermal systems (EGS) reservoir 
conditions. In our previous work (Kalinina et al 2011, 
Kalinina et al 2012a, b) we conducted a number of re- 
servoir simulations assuming homogeneous and hetero- 
geneous reservoir conditions. This work demonstrates 
the importance of natural fracture properties on heat 
extraction. 

The FCM approach was selected because it is compu- 
tationally effective compared to the discrete fracture ap- 
proach. This is especially important when the fracture 
network is very complex and the scale of interest is very 
large, as it is in the case of EGS. A number of different 
techniques were proposed to translate individual frac- 
ture properties into a continuum model (Botros et al. 
2008, McKenna and Reeves, 2006, Reeves et al. 2008). 
The method that we used is an extension of the method 
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described in McKenna and Reeves (2006) and (Reeves 
et al. 2008). The method is based on mapping the fracture 
properties into a continuum model regular grid. The per- 
meability tensor calculated as a part of this mapping 
allows for preserving the properties of fracture net- 
work. A similar approach applied to a dual porosity 
continuum model is proposed in Hao et al. (2013). Hao 
et al. (2013) also demonstrated the ability of the FCM to 
accurately simulate heat transport in fractured systems by 
comparing the modeling results obtained with FCM and a 
discrete fracture model. Our goal was to modify the FCM 
approach to better represent the conditions typical for 
EGS, such as continuous fracture features and high vari- 
ability of fracture density with depth. 

To describe the common EGS conditions, we compiled 
data from the literature on natural fractures in granitic 
rocks (Kalinina et al. 2012a, b). These data suggested that a 
typical EGS reservoir may have 2 to 4 sets of sub-vertical 
fractures with different dips and orientations. Each set may 
have different fracture spacing and aperture. The fracture 
density may also change with depth, and repeating intervals 
with high and low fracture densities are very common. To 
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incorporate these data in the reservoir simulations, we ex- 
tended the previous FCM designed for vertical fractures 
(McKenna and Reeves 2006; Kalinina et al 2012a, b). The 
new FCM is described in (Kalinina et al 2012a, b) and in- 
corporates fully three-dimensional representations of aniso- 
tropic permeability, multiple independent fracture sets, and 
arbitrary fracture dips and orientations. 

The goal of this work is to demonstrate the capabilities 
of the FCM to generate different types of fractures, to 
evaluate the impacts of these different fracture types on 
heat extraction, and to compare alternative heat extrac- 
tion schemes. 

Method 

The FCM maps the permeability of discrete fractures 
onto a regular grid using a continuum approach. This 
approach defines fracture sets using the strike, dip, aper- 
ture, and spacing of parallel fractures in the system. The 
purpose of the FCM is to generate multiple permeability 
fields using field observations of fracture set strike, dip, 
aperture, and spacing. Both, natural fractures or frac- 
tures created with hydro-fracturing can be represented. 

Parallel plate flow methods, originally presented by Snow 
(1968, 1969), were extended to include multiple fracture 
sets at arbitrary fracture orientation following the method 
developed by Chen et al. (1999). 

The continuum approach computes permeability ten- 
sors for each grid cell in the model domain. For one 
fracture set, the permeabflity tensor is defined in Eq. (1) 
(Chen et al. 1999). 
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where kij is the permeability tensor in the / = x,y,z and ; = 
x,y,z directions, b is fracture aperture, d is fracture spacing, 
and n 1^2,3 is the unit normal to the fracture plane in the x, 
y, and z direction, respectively. 

The unit normal components to the fracture plane are 
defined in Eq. (2). 
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where N is the number of fracture sets and is defined 
by Eq. (1). The summation assumes that the total poros- 
ity within a grid-cell changes very little, which is valid 
for most EGS conditions. Note that only /<^^, kyy, and kzz 
components of the permeability tensor are used in the 
heat transport model. 

The kx^T, kyy^ and kzJ^ values are calculated for each 
grid block based on the fracture aperture, spacing, strike, 
and dip probability distributions that are defined for the 
fracture set m. Four probability distributions are currently 
incorporated: normal, log-normal, exponential, and power 
law. However, any other distribution can be easily added. 

SGSIM (Deutsch and Journel, 1998) is used to gene- 
rate a multiGaussian field of spatially correlated num- 
bers. At each grid point, the simulated value can be 
considered a Gaussian deviate transformed to a uniform 
deviate. At each grid point, the uniform deviate is then 
used to draw from a fracture parameter distribution (i.e., 
aperture, spacing, strike, and/or dip). This is essentially 
the probability field approach to stochastic simulation. 
Here, there are no conditioning data, so all distributions for 
the chosen fracture parameter are identical across the do- 
main and the spatial correlation in the simulated uniform 
deviates is translated to the fracture parameter. The orien- 
tation and magnitude of the spatial correlation can be spe- 
cified to create continuity in similarly valued parameters 
in certain directions. Individual fracture sets can utilize 
both correlated/uncorrelated properties (e.g., correlated 
spacing but uncorrected apertures). 

As noted in literature search, fracture properties often 
change with depth. Different fracture parameters can be 
assigned to different model regions to modify fracture 
density with depth. By adjusting the fracture aperture 
and spacing with depth, zones of high and low fracture 
density can be created. 

The resulting /c^^, kyy, and kzz values for each grid cell 
calculated using Eq (1) through Eq. (3) are written in the 
format required for the heat transport simulations with 
the Finite Element Heat and Mass Transfer code FEHM 
(Zyvoloski et al. 1997). 

Applications 

Our goal is to consider FCM applications to common 
EGS conditions and to examine heat extraction as a 
function of these conditions. These conditions include 
the following: 



where a = 90° - dip and co = strike - 90°. 

In the case of multiple fracture sets, the permeabflity 
tensor can be computed by summing the permeabflity 
tensors for individual fracture sets as foflows. 



i. Continuous fractured features. 

ii. Fractures with different spatial orientation 

iii. Multiple fracture sets. 

iv. Depth intervals with different fracture density. 
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The discussion of these applications is provided below. 
The permeability fields were generated for a modeling do- 
main measuring 300 m by 300 m in the horizontal direc- 
tions and 200 m in the vertical direction using a constant 
orthogonal grid with a grid cell size of 5 m by 5 m by 5 m. 

Continuous fractured features 

Natural fractures in granite often form multiple sets of par- 
allel fractures that extend to significant distances. Fractures 
created with hydro-fracturing also may extent to a few 
tenths of meter. Consequently, the FCM needs to have a 
capability to generate continuous fractures. Our first goal 
was to define the fracture parameters or the combination 
of fracture parameters that allows for generating continu- 
ous fracture sets. 

As discussed above, the spatial correlation in fracture 
orientation can be introduced via spatially correlated ran- 
dom numbers, which can then be used to define any com- 
bination of fracture parameters. We first used the spatial 
correlation for one of the fracture parameters while keeping 
the other three parameters uncorrected. We then consid- 
ered different combinations of the spatially correlated frac- 
ture parameters. 

Based on our experiments we concluded that fracture 
spacing is the major parameter controlling the spatial 
continuity of the fractures. The importance of the frac- 
ture aperture and spacing is clear from Eq. (1). However, 
using the spatial correlation for the aperture alone did 
not allow for generating continuous features. 

The following approach to generating continuous frac- 
tured features is recommended. To generate fractures that 



extend through multiple grid cells, fracture spacing is cor- 
related in the direction of strike and dip. The strike, dip, 
and aperture can either use the same correlated field as the 
spacing or a random number generator to assign the value 
for each cell. The probability distribution function for strike 
and dip should be aligned with the spatial correlation for 
spacing. For example, if the spatial correlation for spacing 
has a strike of 45° and a dip of 0°, the mean value for strike 
should be 45° and the mean value for dip should be 0°. In 
this way, fractures extend through multiple grid cells with 
similar properties. The correlation length used for spacing 
controls the length of that fracture set. 

An example is shown in Figures 1 and 2 for one fracture 
set with a mean strike of 135° and a mean dip of 90°. A nor- 
mal distribution with a mean of 3.35x10'^ m (0.335 mm 
and a standard deviation of 3.35x10'^ m (0.0335 mm) was 
used for the fracture aperture. An exponential distribution 
with a mean of 1.0 m, minimum of 0.001 m and maximum 
of 10 m was used for the fracture spacing. The permeability 
field shown in Figure 1 was obtained using spatial correl- 
ation in generating the grid cell fracture apertures, while in 
Figure 2 spatial correlation was used in generating the grid 
cell fracture spacing. Figures 1 and 2 demonstrate that the 
spatial correlation in fracture aperture does not produce con- 
tinuous features while spatial correlation in fracture spacing 
does. Note that in this example the search radius for SGSIM 
was defined to be half of the modeling domain size (150 m). 

Fractures with different spatial orientation 

In our previous work (Kalinina et al. 2012a, b), we demon- 
strated the importance of vertical anisotropy in permeability 
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Figure 1 Permeability field generated using spatial correlation in fracture aperture. Permeability is in 
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Figure 2 Permeability field generated using spatial correlation in fracture spacing. Permeability is in m^. 



for heat extraction using a homogeneous anisotropic reser- 
voir as an example. We showed that heat extraction can be 
significantly improved in the case when horizontal reservoir 
permeability is a few orders of magnitude greater than the 
vertical permeability. We also concluded that vertical frac- 
tures (vertical anisotropy is equal to 1) are the worst case for 
heat extraction because they allow cold water to sink to the 
bottom of the reservoir. 

In this work we examine fractures with different dips, 
evaluate the average vertical anisotropy as a function of 
fracture dip, and demonstrate the impacts of spatial orien- 
tation on heat extraction. 



The average ratio of horizontal and vertical permeability 
(kxx/kzz or kyy/kzz) as a function of mean fracture dip can 
be calculated from Eq. (1) and Eq. (2) using the mean frac- 
ture parameters. Note that in the case when the strike is ei- 
ther 45° or 135°, kxx = kyy and kxx/kzz = kyy/kzz- Figure 3 
demonstrates kxx/kzz as a function of fracture dip for the 
case with the mean fracture strike of 135°, mean aperture 
of 3.35x10'^ m, and mean spacing 1.0 m. This functional 
relationship is very similar to the other cases with different 
mean fracture spacing and aperture values. 

The vertical anisotropy is 2 orders of magnitude or greater 
(Figure 3) only in the case of sub-horizontal fractures with a 
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Figure 3 Vertical anisotropy as a function of fracture dip. 
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Figure 4 Permeability field for one set of horizontal fractures. Permeability is in 
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dip of less than 6°. In the case of sub-vertical fractures 
(dip greater than 60°), the vertical permeability is greater 
than the horizontal permeability. Note that natural frac- 
tures in deep granite rocks are mostly sub-vertical, which 
leads to insignificant anisotropy in vertical direction. 

Heat extraction for different fracture settings 

To investigate in impact of various fracture properties 
on heat extraction simulations, the following cases were 
considered: 



i. One set of close to vertical fracture (strike of 135°, 
dip 85°) 

ii. One set of close to horizontal fractures (strike of 
135°, dip 2°) 

iii. One set of fractures with dip of 10° (strike of 
135°) 

iv. Two sets of close to horizontal nearly orthogonal 
fractures (strike of 130° and 45°, dip 2°) 

V. Different fracture sets within the different depth 
intervals 
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Figure 5 Permeability field for two sets of horizontal fractures. Permeability is in m^. 
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Figure 6 Permeability field for one set of horizontal fractures with different aperture and spacing within different depth intervals. 

Permeability is in m^. 



The close to horizontal fractures are not common in gra- 
nite rocks at large depth. However, these fractures might be 
created with hydro-fracturing. This was one purpose for 
considering the close to horizontal fractures in the heat ex- 
traction simulations. Another purpose was to evaluate the 
effects of vertical anisotropy, which is at its maximum value 
for horizontal fractures. 



The permeability fields for these examples were gener- 
ated using FCM as defined by Eq. 1 through 3. Figure 2 
shows the permeability field generated for the vertical 
fracture case. Figures 4, 5, and 6 show the permeability 
fields for the other cases, except the case with a dip of 
10°, which is similar to the case with horizontal frac- 
tures. These figures demonstrate that FCM is capable of 
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Figure 8 Average temperature in the production well for different fracture settings. 



generating continuous fractured features for a range of 
conditions. 

The thermo-hydrology simulations in a fractured res- 
ervoir were conducted using the numerical computer 
code FEHM (Zyvoloski et al. 1997). FEHM is a control 
volume finite element code for simulating subsurface 
non-isothermal multi-phase multi-fluid heat and mass 
transfer developed at Los Alamos National Laboratory 
(LANL). The code has been extensively used for model- 
ing geothermal systems and groundwater flow and trans- 
port, including large scale projects, such as performance 
assessment of Yucca Mountain Project and Environmen- 
tal Remediation of the Nevada Test Site. The governing 
equations for the conservation of mass and energy im- 
plemented in FEHM are provided below. More details 
can be found in (Zyvoloski et al. 1997). 

Conservation of mass for water is described as: 

dA 

-^+V.fm+q„ = 0 (4) 

where the mass per unit volume, A^, is given by: 

Am = (p{Sy p^{l-rj^) +Sipi{l-rji)) (5) 

and the mass flux, fm is given by: 

fm = {l-%)p^Vv + {^-rji)vipi 

where, (|) is the porosity of the matrix, S is saturation, p 
is density, rj is the concentration of the noncondensible 
gas and is expressed as a fraction of the total mass, v is 
velocity, and the subscripts v and / indicate quantities 
for the vapor phase and the liquid phase, respectively. 
Source and sink terms are represented by the term q^. 
Conservation of fluid- rock energy is described as: 



-^+V-f, + q, = 0 (6) 
where the energy per unit volume, A^, is given by: 

Ae = (l-(p)p^Ur + Cp {SypP^Uy + SipiUl), 

Ur = CprT (J) 

and the energy flux, is given by: 

/; =P^hyVy^VlPlhl-KT (8) 

The subscript r refers to the rock matrix, Uy, Uy, and Ui 
are specific internal energies, c^^ is the specific heat, h^, 
and h-i are specific enthalpies, K is an effective thermal 
conductivity, T is the temperature, and qe is the energy 
contributed from sources and sinks. 

It is assumed that Darcys Law applies to the move- 
ment of each phase and the velocity can be expressed as: 



Vy=^{yPy-p^-g) (9a) 

vi=^{vPi-Pig) (9b) 

Py=Pl+Pcap (9c) 



where k is the permeability, Ry, and Rj^ are the relative per- 
meabilities, py and pi are viscosities, Py and Pi are the 
phase pressures, Pcap is the capillary pressure and , g rep- 
resents the acceleration due to gravity. The equations are 
shown for an isotropic medium, though this restriction 
does not exist in the computer code. 

We used the FEHM capability to simulate heat and 
mass transfer with pressure and temperature dependent 
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properties required for modeling a geothermal system. 
FEHM incorporates constitutive relationships that describe 
pressure and temperature dependent fluid and gas proper- 
ties, relative permeabilities and capillary pressures, stress 
dependencies, and reactive and sorbing solutes. 

The considered model (Figure 7) represents H of the 
reservoir. This was done to take advantage of a symmetry 
introduced by a 5 -point injection scheme in which there is 
one injection well in the middle and 4 production wells in 
the corners. The model is 300 m by 300 m in the horizontal 



directions and 200 m in the vertical direction. The grid 
block size is 5 m by 5 m by 5 m. The total production 
rate was 30 kg/s. The injection well was located in the cor- 
ner of the modeling domain at = 0, j = 0 m. The produc- 
tion well was located 417 m away from the injection well 
along the model diagonal {x = 300 m, 3/ = 300 m). The 
fracture trend was aligned with the model diagonal con- 
necting the injection and the production well. The ini- 
tial reservoir temperature was 225°C and the injection 
temperature was 80°C. 
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Figure 1 1 Temperature Distribution after 30 Years of Injection for Reservoir with Different Fracture Density in Different Depth Intervals. 



The pre-injection boundary conditions are hydrostatic 
pressures corresponding to the depth to the reservoir top 
(4 km) and constant temperature on the vertical sides of 
the model, constant pressure and temperature on the top 
of the model, and zero flux on the bottom of the model 
Zero flux conditions are specifled across all the boundaries 
during the injection. Injection is implemented using a fixed 
injection mass-flow rate. The production well is imple- 
mented by specifying a fixed bottom-hole pressure in the 
production well. 

Figure 8 shows the results of the heat simulations for the 
different fracture settings in terms of average temperature 
in the production well as a function of time from the begin- 
ning of injection. 

Figure 8 demonstrates that heat extraction in the case 
with fractures dipping 10° (kxx/kzz = 33) is similar to the 
case with the vertical fractures. The temperatures at the 
end of the injection are the same in both cases. This result 
is consistent with our previous work in which we consid- 
ered a homogeneous anisotropic reservoir (Kalinina et al. 
2012a, b). 

The heat extraction in the case of horizontal frac- 
tures (kxx/kzz =3,119) is significantly better than in the 
case of vertical and sub-horizontal (dip 10°) fractures. 
The temperature begins to decrease 5 years later and 
is 10 °C higher by the end of injection. There is a small 
difference in performance (temperature histories) be- 
tween the case with one set of horizontal fractures and 
two sets of horizontal fractures. 

The temperature distribution at the end of the si- 
mulation is shown in Figure 9 for the vertical fracture 



case and in Figure 10 for the horizontal fracture case. 
These figures show that a larger reservoir volume is 
in the heat exchange in the horizontal fracture case. 
This leads to slower cooling and higher production 
temperatures. 

The heat extraction is also improved in the case when 
there are intervals with high density and low density hori- 
zontal fractures compared to the case with horizontal frac- 
tures with constant density. The production temperature 
after 10 years is 3° to 6° higher. However, the higher 
temperature were only obtained in the case when the verti- 
cal permeability in the low fracture density intervals is 2 or- 
ders of magnitude lower than in the high fracture density 
intervals. Otherwise, the production temperatures are simi- 
lar in both cases. The intervals with low fracture density 
and low vertical permeability create additional vertical an- 
isotropy. As a result, the temperature front is spread over a 
larger reservoir volume, which improves heat exchange 
and lead to slower cooling and higher production tempera- 
tures. The temperature distribution at the end of injection 
for this case is shown in Figure 11. 

Alternative heat extraction schemes 

Our simulations for the different fracture settings dem- 
onstrated that the best heat extraction conditions are in 
the case of horizontal fractures. Even gently dipping 
fractures impact the heat extraction in a negative way 
and make it simflar to the worst case of sub-vertical 
fractures. Because the natural fractures in granite rock 
are sub-vertical, the extraction of heat with vertical wells 
will likely be limited. 
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One solution to this problem might be to use a differ- 
ent heat extraction scheme. An example of an alternative 
scheme is considered below. 

The proposed alternative scheme uses horizontal instead 
of vertical wells. Both schemes are shown in Figure 12. To 
make these schemes comparable, well length, well separ- 
ation and total production rate need to be the same in 
both cases. Two options were considered for the horizon- 
tal well scheme: the injection well below the production 
well, and the injection well above the production well. 

The vertical well case uses one injector and 2 producers, 
which is a 3-point injection scheme. The horizontal well 
case uses two injectors and two producers, which is a 4- 
point injection scheme. In both cases the total production 
is 15 kg/s, the well length is 200 m and the well separation 
is 150 m. The H modeling domain is used to simulate heat 



extraction as shown in Figure 12. The only difference be- 
tween the two cases is in the well orientation. The simula- 
tions were done for a homogeneous isotropic reservoir 
with permeability 1x10'^^ m^ and for a reservoir with one 
set of vertical fractures similar to one shown in Figure 1, 
but with the fracture trending along the y axis. 

Figure 13 shows the results of the heat simulations for 
the vertical and horizontal well configurations in terms 
of average temperature in the production well as a func- 
tion of time from the beginning of injection. 

Figure 13 indicates that the heat extraction for the hori- 
zontal wells with the injector at the bottom is significantly 
better than the heat extraction with the vertical wells. As 
was expected, placing the injector at the top in the hori- 
zontal well scheme negatively impacts the heat extraction 
as compared to the case with the injector at the bottom. 
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Figure 13 Average Temperature in the Production Well for Alternative Heat Extraction Schemes. Heterogeneous case refers to a reservoir 
with the vertical fractures. 



The heat extraction in the case with the injector at the top 
is worse than in the case with the vertical wells. 

The performance of the reservoir with the vertical 
fractures is slightly better than the performance of a 
homogeneous reservoir in the case of vertical wells. This 
is consistent with our previous conclusions (Kalinina 
et al. 2012a, b). 



The opposite effect was observed in the case of hori- 
zontal wells. The temperature distribution in reservoir 
is shown in Figure 14 for the homogeneous reservoir 
and in Figure 15 for the reservoir with one set of verti- 
cal fractures. The vertical fractures resulted in perco- 
lating temperature front in the vertical direction. This 
negatively affected the heat extraction in the vertical 
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fracture case compared to the homogeneous reservoir 
case. 

Conclusions 

The fractured continuum model (FCM) provides a tool for 
mapping the permeability of discrete fractures (either 
natural or created with hydro- fracturing) onto a regular 
modeling grid using a continuum approach. The fracture 
characteristics obtained in the field, such as number of frac- 
ture sets, fracture orientation, spacing and aperture, can be 
used to generate multiple reservoir representations. The 
FCM can be extended to represent continuous fracture fea- 
tures and fracture intervals with different fracture densities. 

The FCM was used to generate permeability fields repre- 
senting different fracture properties. A 5-point injection 
scheme was considered for each representation to evaluate 
the impacts of fracture properties on heat extraction. The 
following conclusions were made based on the heat trans- 
port simulations. 

■ The vertical anisotropy is the major factor affecting heat 
extraction. The heat extraction in the case with fractures 
dipping 10° (kxx/kzz = 33) is only slightly better than in the 
case with the vertical fractures. 

The heat extraction in the case of horizontal fractures (k^x/ 
kzz = 3,119) is significantly better than in the case of vertical 
and sub-horizontal (dip 10°) fractures. The improved heat 
performance in the horizontal fracture case is related to a 
larger reservoir volume involved in the heat exchange. 

■ The average vertical anisotropy is a function of fracture 
dip. The vertical anisotropy is significant (2 orders of 



magnitude or greater) only in the case of sub-horizontal 
fractures with a dip less than 6°. In the case of sub- 
vertical fractures (dip greater than 60°), the vertical per- 
meability is greater than the horizontal permeability. 

■ The heat extraction is further enhanced in the case 
when there are intervals of high fracture density and low 
fracture density as compared to the case with constant 
fracture density with depth. However, the performance is 
only noticeable better when the vertical permeability in 
the low density intervals is 2 orders of magnitude lower 
than in high density intervals. This suggests that creating 
fracture zones with intact rock between them may be a 
good way to design an efficient reservoir. 

■ The heat extraction with the horizontal wells with the 
injector at the bottom is significantly better than the heat 
extraction with the vertical wells. In the case with 
horizontal wells, the temperature remains constant for 
the first 15 years of injection and then starts dropping at 
a rate of 0.7 degree C/year from 15 to 20 years and at the 
rate of 1.9 degree C/year during the remaining injection 
time. In the case with vertical wells, the temperature 
drops nearly linearly from the very beginning of injection 
at a rate of 2.2 degree C/year. 

Competing Interests 

The authors declare that they have no competir^g interests. 
Authors' contributions 

EK and TH developed the heat and mass transfer model and performed the 
simulations. KK and SM developed an extension to the fractured continuum 
approach. EK, TH, and TL developed the simulation cases and analyzed the 
results. All authors read and approved the final manuscript. 



Kalinina et a I. SpringerPlus 2014, 3:1 10 
http://www.springerplus.conn/content/3/1 /1 1 0 



Page 13 of 13 



Acknowledgments 

Sandia National Laboratories is a multi-program laboratory managed and 
operated by Sandia Corporation, a wholly owned subsidiary of Lockheed 
Martin Corporation, for the U.S. Department of Energy's National Nuclear 
Security Administration under contract DE-AC04-94AL85000. 

Author details 

^Sandia National Laboratories, P.O. Box 5800, Albuquerque, NM 87185, USA. 
^IBM Research, Smarter Cities Technology Centre, BIdg. 3, Damastown 
Industrial Estate, Mulhuddart, Dublin 15, Ireland. 

Received: 21 January 2014 Accepted: 14 February 2014 
Published: 24 February 2014 

References 

Botros FE, Hassan AE, Reeves DM, Pohll G (2008) On mapping fracture networks 
onto continuum. Water Resour Res vol 44:W08435 

Chen M, Bai M, Roegiers J-C (1999) Permeability tensors of anisotropic fracture 
networks. Math Geol 31:4 

Deutsch CV, Journel AG (1998) Gslib Geostatistical Software Library And User's 
Guide, 2nd edn. Oxford Univ Press, New York 

Hao Y, Fu P, Carrigan CR (2013) Application Of A Dual-Continuum Model For 
Simulation Of Fluid Flow And Heat Transfer In Fractured Geothermal 
Reservoirs. Proceedings, 38th Workshop On Geothermal Reservoir Engineering, 
vol SGP-TR-198. Stanford University, Stanford, California, pp 462-469 

Kalinina E, Hadgu T, McKenna SA, Lowry TS (201 1) Bridging the gap between 
complex numerical modeling and rapid scenario assessment: a 
dimensionless parameter approach for enhanced geothermal systems. 
Trans Geoth Resour Counc 35:407-412 

Kalinina E, McKenna SA, Hadgu T Lowry TS (2012a) Analysis Of The Effects Of 
Heterogeneity On Heat Extraction In An Egs Represented With The 
Continuum Fracture Model. Proceedings, 37-th Workshop on Geothermal 
Reservoir Engineering. Stanford University, Stanford, California, pp 436-445, 
January 30 - February 1, 2012, SGP-TR-194 

Kalinina E, McKenna SA, Klise K, Hadgu T, Lowry TS (2012b) Incorporating 
complex three-dimensional fracture network into geothermal reservoir 
simulations. Trans Geoth Resour Counc 36:493-498 

McKenna SA, Reeves PC (2006) Fractured Continuum Approach To Stochastic 
Permeability Modeling. In: Coburn TC, Yarus JM, Chambers RL (eds) 
Stochastic Modeling and Geostatistics: Principles, Methods, and Case Studies, 
Volume II. American Association of Petroleum Geologists, Tulsa, Oklahoma, 
pp 173-186 

Reeves DM, Benson DA, Meerschaert MM (2008) Transport of conservative solutes 
in simulated fracture networks: 1. Synthetic data generation. Water Resour 
Res vol 44:W05404 

Snow DT (1968) Rock fracture spacings, openings, and porosities. J Soil Mech 

Foundations Div, Proc Am Soc Civ Eng 94:73-91 
Snow DT (1969) Anisotropic permeability of fractured media. Water Resour Res 

5:1273-1289 

Zyvoloski GA, Robinson PA, Dash ZV,Trease LL (1997) Summary Of The Models 
And Methods For The Fehm Application - A Finite-Element Heat- And 
Mass-Transfer Code. Los Alamos National Laboratory, LA-13307-MS, 
Los Alamos, New Mexico 



doi:1 0.1 186/21 93-1 801 -3-1 10 

Cite this article as: Kalinina et al.: Applications of fractured continuum 
model to enhanced geothermal system heat extraction problems. 

SpringerPlus 2014 3:110. 



f ^ 

Submit your manuscript to a SpringerOpen^ 
journal and benefit from: 

► Convenient online submission 

► Rigorous peer review 

► Immediate publication on acceptance 

► Open access: articles freely available online 

► High visibility within the field 

► Retaining the copyright to your article 



Submit your next manuscript at ► springeropen.com 



J 



